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Abstract 

Background: Information about the genotypic characteristic of Coxiella burnetii from Hungary is lacking. The aim 
of this study is to describe the genetic diversity of C. burnetii in Hungary and compare genotypes with those found 
elsewhere. A total of 12 samples: (cattle, n = 6, sheep, n = 5 and human, n = 1) collected from across Hungary were 
studied by a 10-loci multispacer sequence typing (MST) and 6-loci multiple-locus variable-number of tandem repeat 
analysis (MLVA). Phylogenetic relationships among MST genotypes show how these Hungarian samples are related 
to others collected around the world. 

Results: Three MST genotypes were identified: sequence type (ST) 20 has also been identified in ruminants from 
other European countries and the USA, ST28 was previously identified in Kazakhstan, and the proposed ST37 is 
novel. All MST genotypes yielded different MLVA genotypes and three different MLVA genotypes were identified 
within ST20 samples alone. Two novel MLVA types 0-9-5-5-6-2 (AG) and 0-8-4-5-6-2 (AF) (Ms23-Ms24-Ms27-Ms28- 
Ms33-Ms34) were defined in the ovine materials correlated with ST28 and ST37. Samples from different parts of 
the phylogenetic tree were associated with different hosts, suggesting host-specific adaptations. 

Conclusions: Even with the limited number of samples analysed, this study revealed high genetic diversity 
among C. burnetii in Hungary. Understanding the background genetic diversity will be essential in identifying and 
controlling outbreaks. 
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Background through the consumption of raw milk and milk products 

The aetiological agent of Q (Query) fever is Coxiella [2]. In humans, Q fever is typically an acute febrile illness 
burnetii, a Gram-negative, obligate intracellular, zoonotic with nonspecific clinical signs, such as atypical pneumonia 
bacterium which is widely distributed throughout the and hepatitis in roughly 40% of the cases while the 60% of 
world. Domestic ruminants (cattle, sheep and goats) the individuals remain asymptomatic after infection [3]. A 
have the main role in the disease cycle, but C. burnetii small percentage (~5%) of infected people may develop 
has been isolated from other domestic and wild animals chronic infection with life-threatening valvular endocar- 
too [1]. Domestic ruminants are reservoirs of the agent ditis [4,5]. 

with usually subclinical manifestation of the disease, but Q fever is a notifiable disease in Hungary. Antibodies 
may suffer from reproductive disorders, and abortion against C. burnetii were first detected in the sera of abattoir 
and stillbirth can occur. The bacteria are transmitted workers in 1950 [6], and infections were first diagnosed in 
mainly by inhalation. Additionally, infections may occur 1956 in dairy and sheep farms [7]. According to recently 

published data, seroprevalence among cattle and sheep 
in Hungary were 38.0% and 6.0% with enzyme-linked 
immunosorbent assay respectively, which correspond 
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reported acute human infections in Hungary ranged be- 
tween 36 and 68 during the past five years. Major outbreaks 
were last registered in the period of 1976-1980 and in 2013 
in Hungary (unpublished data). 

Genetic characterisation of C. burnetii is required for 
epidemiological investigations in Q fever outbreaks and 
for surveillance purposes. Several typing systems exist, 
including pulsed-field gel electrophoresis [10], sequence 
analysis or restriction fragment length polymorphism of 
single genes (16S ribosomal RNA, icd, coml or rawcZ) 
[11-15]. The main problems with these methods are that 
their discriminatory power is poor and their reproduci- 
bility and transferability is not always straightforward 
[16]. Recently, two PCR-based typing methods have 
come into view: multispacer sequence typing (MST) [17] 
and multi-locus variable-number tandem repeat analysis 
(MLVA) [18,19]. Both methods have high discriminatory 
power and are reproducible, allowing interlaboratory com- 
parisons without the need for isolation and cultivation of 
the organism under biosafety level three conditions. 

Multiple sequence types have been identified from the 
Central European region in ticks and human samples: ST2 
from Ukraine, ST16, ST18 from Romania and Slovakia, 
ST23 from Czech Republic and Slovakia and ST22, ST29 
from Slovakia. While ST27 and ST32 were detected in 
goat samples in Austria [17,20]. MLVA types L, P and par- 
tial MLVA profiles were identified from cow milk samples 
originating from Slovakia, Austria and Croatia [21]. 

The aim of this study was to identify C. burnetii MST 
and MLVA genotypes occurring in Hungary and to com- 
pare them with types from foreign countries. 



Methods 

This study included five ovine and five cattle abortion 
samples (cotyledons) and a cattle milk sample collected 
from different parts of Hungary (Table 1, Figure 1). The 
blood sample of a 72 years-old man suffering from acute 
Q fever (1:128 IgGII titre with micro-immunofluorescence 
test, Focus Diagnostics, Cypress, CA) was also character- 
ized. The samples were not related to any outbreak, they 
were collected through routine diagnostic examinations. 
All tests were performed in accordance with all applicable 
institutional and national guidelines and regulations, 
approved by the ethics committees of the Institute for 
Veterinary Medical Research (animal) and the National 
Center for Epidemiology (human) and with the consent of 
the patient. 

Extraction of DNA from the samples was performed 
with the QIAamp Tissue and Blood kit (Qiagen GmbH, 
Hilden, Germany) according to the manufacturers instruc- 
tions. To quantify the approximate bacterial load and to 
determine the cycle threshold (Ct) values, all samples were 
analysed by TaqMan based real-time polymerase chain 
reaction (PCR) assay targeting the multicopy IS1111 inser- 
tion element of the C. burnetii genome [22] . 

MST consisted of amplification followed by sequencing 
of ten different spacer regions of the C. burnetii genome: 
Cox2, 5, 18, 20, 22, 37, 51, 56, 57 and 61 [17]. Primer 
sequences and reaction conditions were described earlier 
[17]. The compiled sequences of the PCR products were 
analysed with the BioEdit Sequence Alignment Editor 
7.1.11 [24]. Sequence types (STs) were determined using 
the MST database and previous publications [17,20,25]. 



Table 1 MLVA and MST genotypes of C. burnetii specimens detected from cattle, sheep and a human clinical sample in 
Hungary 



Sample ID 


Source 


Town of origin 


Year 


Ct 


Ms23 


Ms24 


MLVA-6 
Ms27 Ms28 


Ms33 


Ms34 


MLVA type 


MST type 


Coxi-14 


cow milk 


Martonvasar 


2010 


32.90 


6 


13 


2 


7 


9 


9 


1 


ST20 


VS41 


cow cotyledon 


Nagyecsed 


201 


2 


2778 


6 


13 


2 


/ 


9 


9 




ST20 


VS76 


cow cotyledon 


Szombathely 


201 


2 


14.35 


6 


13 


2 


/ 


9 


10 


J 


ST20 


VS79 


cow cotyledon 


Nyirbator 


201 


2 


30.34 


6 


13 


2 


7 


9 


10 


J 


ST20 


VS93 


cow cotyledon 


Berkesd 


201 


2 


15.39 


6 


13 


2 


/ 


9 


11 


M 


ST20 


VS103 


cow cotyledon 


Kapuvar 


201 


2 


16.89 


6 


13 


2 


7 


9 


9 




ST20 


VS16 


sheep cotyledon 


Balkany 


201 


2 


8.85 


0 


8 


4 


5 


6 


2 


AF 


ST37 


VS27 


sheep cotyledon 


Jaszfelsoszentgyorgy 


201 


2 


29.83 


0 


8 


4 


5 


6 


2 


AF 


ST37 


VS38 


sheep cotyledon 


Beremend 


201 


2 


12.79 


0 


8 


4 


5 


6 


2 


AF 


ST37 


VS42 


sheep cotyledon 


Kunmadaras 


201 


2 


30.72 


0 


8 


4 


5 


6 


2 


AF 


ST37 


VS108 


sheep cotyledon 


Biharnagybajom 


201 


3 


8.73 


0 


9 


5 


5 


6 


2 


AG 


ST28 


4756 


human blood 


Hortobagy 


201 


1 


28.31 


5 


9 


5 


ND 


6 


2 


partial 


ND 


Cburnetii RSA493* 


tick 


USA 


1935 




9 


27 


4 


6 


9 


5 




ST16 



Ct: cycle threshold value of real-time PCR assay targeting the multicopy IS 7 7 1 7 insertion element of the C. burnetii genome [22] to quantify the approximate 
bacterial load, 0: no amplification, ND: not determined probably because of insufficient DNA content, *: Genotype of RSA493 strain (Nine Mile, GenBank accession 
number: AE016828) was determined by in silico analysis. 
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Figure 1 Map of Hungary showing the geographical distribution of C. burnetii from cattle (red circles), sheep (blue triangles) and human 
(green star). MST genotypes and MLVA types of samples are also presented. ND: not determined. (The blank map was downloaded from an open 
source [23]). 



Unfortunately, the allele codes for ST35 and ST36 are 
not publicly available therefore we were not able to 
include them in the phylogenetic analysis [25]. To deter- 
mine the phylogenetic relationships of our samples to 
known MST genotypes [17,20] we used the sequence 
information of our three STs and constructed a phylo- 
genetic tree with the 112 polymorphisms and methods 
described in Hornstra et al. [26]. As samples VS16/27/ 
38/42 (sheep abortion) appeared to have a new ST and 
homoplasies can alter topologies of trees, we also con- 
structed a sub-tree to verify its placement using the 
same methods and 112 polymorphisms as above but only 
with STs: 8, 9, 10, 27, 28, 31, and the two sheep abortion 
STs from this study. 

MLVA was performed by single PCRs targeting six 
variable microsatellite markers as performed in other 
recent studies [18,27-30]. Ms27, Ms28 and Ms34 contain 
repeat units of six base pairs and Ms23, Ms24 and Ms33 
contain repeat units of seven base pairs. The 3 '-end-la- 
belled forward and reverse primer sequences, and PCR 
conditions were applied as described before [27-30]. The 
amplification products were run on an ABI 3100 Genetic 
Analyser and electropherograms were evaluated with the 
Peak Scanner Software 2.0 (Applied Biosystems Inc., 
Foster City, CA). DNA of the Nine Mile strain (RSA 
493, Coxevac, Ceva Inc., Budapest, Hungary) was used 
as a reference. The repeat numbers of each marker were 
determined by extrapolation using the obtained length of 
the sample fragments relative to the obtained fragment- 
length of the reference strain. A modified coding conven- 
tion was used for Ms33: nine imperfect repeat units were 



counted in the Nine Mile strain (G. Vergnaud, personal 
communication). Novel MLVA types were determined if 
the described allele combinations had not been previously 
described in any publications applying the same MLVA 
method [21,27-29]. 

Results 

All 12 samples were positive with the Coxiella specific 
real-time PCR with Ct values below 33 and all but one 
yielded complete MST and MLVA genotypes (Table 1). 
MST alignment and MLVA profiles were deposited in 
Dryad ([31], doi:10.5061/dryad.rc8q4) and in the Cox- 
iella MLVA database ([32], accession numbers: HungOO- 
Hungll). 

The MST analysis of the six cattie and five ovine sam- 
ples revealed a novel (proposed ST37) and two previously 
described (ST20, ST28) STs (Table 1, Figure 2). The allele 
codes of the newly identified genotype ST37 are 5-4-9-5- 
8-5-2-3-4-6 for spacers Cox2-Cox5-Coxl8-Cox20-Cox22- 
Cox37-Cox51-Cox56-Cox57-Cox61, respectively. Phylo- 
genetic analysis according to Honstra et al. [26] resulted 
in 156 equally parsimonious tree topologies, one of which 
is shown in Figure 2, with a homoplasy index (excluding 
uninformative characters) of 0.0769. Sub-tree construction 
using STs 8, 9, 10, 27, 28, 33, and the two sheep abortion 
STs from this study resulted in a single tree with no 
homoplasy whose topology exactly matched the topology 
for this clade in the tree selected in Figure 2. In clonal 
groups of organisms with little genetic diversity the homo- 
plasy index has been shown to be the most appropriate 
measure of phylogenetic accuracy [33,34], therefore the 
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iSTI-human-FR 



ST3-human-FR 
ST2-human-FR, UA, KG, CF 
- ST4-human-FR, ES, PT 
ST30-goat-NA 
ST5-human-FR 
ST6-human, tick-FR, SN 
ST7-human-FR, RU 



ST27-goat-AT 

ST31-human-UZ 



— ST8-human, goat, sheep-FR, ES, US, PT, CA 
ST9-human, pigeon-FR 
I — ST10-human, sheep-FR 
ST28-sheep, cattle, tick, human-KZ 



sheep abortion-VS108-MLVA type AG 



- STDugway 



ST21 -human, cat, dog-FR, CA, US 
ST12-human, goat, vegetation-FR, CH, DE 
ST13-human, cattle, goat, sheep, vegetation-FR, ES, PT 
STll-human-FR 
ST34-FR 

ST33-human, sheep, goat-NL, DE, FR 
ST14-human, sheep-FR, DE 
ST15-human-FR 
i— ST24-cattle-DE 
^ ST32-goat-AT 
ST18-human, goat, sheep, cattle-FR, IT, SK, RO, GR, DE, ES, PL 
ST25-fly, vaccine, human-DE, RU, KZ, IT 
i— ST22-tick-SK 

-\ ST23-tick, cattle, bird, mouse, human, bed bug, sheep-RU, UZ, CZ, SK, KZ, MN 

I ST29-human-IT, SK 

- ST19-human-SN 

, ST20-human, cattle, rodent, goat, soil-FR, DE, US, ES, NL, GB 



cattle abortion, milk-VS41/76/79/93/103. Coxi14-MLVA types l/J/M 



■ ST17-human-FR, GF 
I ST16-human, cattle, sheep, tick-FR, US, RO, CF, SK, DE, JP, PL 
- ST26-human-UZ 

Figure 2 Parsimony tree showing the placement and phylogenetic relationships of the sequence types (ST) from this study with 
known STs [19], rooted according to Pearson et al. [33]. Highlights indicate the placement of the STs of our samples; the red highlight 
denotes the newly proposed ST37. One hundred and twelve polymorphisms from Hornstra et al. [26] were used to construct this phylogeny and 
resulted in 156 equally parsimonious trees, the first of which is depicted above. Tree length is 120 and the homoplasy index (excluding 
uninformative characters) is 0.0769. Isolate origins and sources are given according to previous publications [19,20,25,26,28,29,35-38] using the 
following location codes: Austria (AT), Canada (CA), Central African Republic (CF), Czech Republic (CZ), France (FR), French Guiana (GF), Germany 
(DE), Greece (GR), Italy (IT), Japan (JP), Kazakhstan (KZ), Kyrgyzstan (KG), Mongolia (MN), Namibia (NA), Netherlands (NL), Poland (PL), Portugal (PT), 
Romania (RO), Russian Federation (RU), Senegal (SN), Slovakia (SK), Spain (ES), Switzerland (CH), Ukraine (UA), United Kingdom (GB), United States 
(US) and Uzbekistan (UZ). 



results of this sub-tree provide confidence in the accuracy 
of the placement of this novel ST (proposed ST37) in 
Figure 2. 

Five MLVA genotypes were identified among the 11 
cattle and ovine samples with complete genotypes, out 
of which two novel MLVA types were described in this 
study. Three MLVA types (I, J, M) were detected among 
the six cattle samples, and the two novel genotypes 
(AF and AG) were defined from the five ovine samples. 
In the case of the human sample only a partial genotype 
was obtained despite the relatively high DNA load (Ct 
value was 28.31). The identified allele combination in the 
human sample did not match any previously defined 
genotypes, implying the discovery of a novel MLVA type 
as well. The partial MLVA genotype from the human 
sample showed high similarity to the novel genotype AG 



defined in sheep, as four out of the five microsatellites 
were the same. 

All samples originating from cattle were identified as 
ST20, regardless of their MLVA types. Samples from 
sheep belonged to two MST types: ST28 was detected in 
a sample with AG MLVA genotype, and the proposed 
ST37 was described in samples with AF MLVA type. 
Host specificity of the genotypes was observed as all 
MLVA genotypes from ST20 were described from cattle 
only and the two other genotypes were detected in sheep 
exclusively. 

Discussion 

Molecular characterisation of C. burnetii is a useful tool 
to explore the genotypic diversity in a region and to 
determine relationships between variants of this bacterium. 
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The main advantage of both MST and MLVA typing 
methods is their high discriminatory power and that 
the cultivation of the bacteria is unnecessary for their 
application. 

MST and MLVA genotyping methods revealed geno- 
typic diversity among samples from domestic ruminants 
and human in Hungary. Three C. burnetii STs (ST20, 
ST28 and ST37) are present in Hungarian ruminants. 

Indicating the broad range of ST20 occurrence, other 
samples have been reported from The Netherlands, 
Germany, Spain, in human clinical samples from France 
and in cow's milk, soil, and goat placental material from 
the United States [17,26,29,35,39]. Microvariants of the 
ST20 genotype occur as well and have been found 
elsewhere; MLVA genotypes I, J and M were described 
before from cow milk and milk products from several 
countries in Europe (France, the Netherlands, Portugal, 
Spain and Switzerland) and other parts of the world (e.g. 
Qatar, Saudi Arabia) [21]. ST20 was also detected in an 
outbreak among goats recorded from the United Kingdom 
[36]. In conclusion, although ST20 is most often asso- 
ciated with cattle products, this genotype occasionally 
infects other species as well. 

ST28 samples have been collected from cows, sheep, 
Hyalomma sp. tick and human blood in Kazakhstan 
[17]. ST37 genotype is closely related to ST27-28, 31 
and 8-10 genotypes that have been collected from 
humans in Austria, Kazakhstan, Uzbekistan, France and 
Portugal as well as sheep and goats from Spain, France, 
Canada and the USA (Figure 2) [17,20,28,29,39]. The 
comparison of the novel genotypes (AG and AF) with 
others described in Europe revealed similarities with geno- 
types AA and T, differing at 2 and 3 loci respectively and 
also associated with sheep, goats and human [28,29]. 
Again, that these genotypes are most commonly asso- 
ciated with sheep and goats provides further evidence for 
host preference [39]. We speculate that there may be 
some association between the Hungarian human and 
ovine infections (MLVA genotype AG) given genotype 
similarities and as the geographical distance between the 
two samples was only 50 km. 

No matches were found among the MST and MLVA 
profiles from Hungary and from the dominant genotype 
(ST33) associated with the Q fever outbreak of The 
Netherlands that emerged in 2007 [27,35]. STs described 
earlier from other Central European countries differed 
also from the STs found in Hungary [17,20]. 

Conclusions 

The present study provides information about the geno- 
typic diversity of C. burnetii occurring in Hungary. From 
the examined 12 C. burnetii samples originating from 
cattle, sheep and human three MST types (containing 
one novel profile) and five MLVA types (including two 



novel profiles) have been identified. These genotypes 
cluster with host, suggesting host-specific evolutionary 
adaptations as shown by Pearson et al. [39] as well. The 
additional analyses of samples from a broader range of 
livestock species are needed to confirm this finding. Our 
results provide data for the European surveillance and 
help to increase information about Q fever in this Central 
European region. 
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